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Abstract. Brief review is given on recent numerical research of the ordering of two typical 
models of spin glasses (SGs), the three-dimensional (3D) Ising SG and the 3D Heisenberg SG 
models. Particular attention is paid to the questions of whether there is a thermodynamic 
transition in zero field, what are the associated critical properties, what is the nature of the 
ordered state, particularly of a possible replica-symmetry breaking, and whether there is a 
thermodynamic transition in applied fields. The properties of the two models are contrasted, 
and possible relation to experiments is discussed. 



1. Introduction 

Spin glasses are the type of random magnets in which both ferromagnetic and antiferromagnetic 
interactions coexist and compete, thereby giving rise to the effects of frustration and quenched 
randomness. Spin glasses (SGs) have been regarded as a typical example of complex systems, and 
their ordering properties have been studied quite extensively, either experimentally, analytically 
or numerically [T]. Experimentally, we now have fairly convincing evidence that SG magnets 
exhibit an equilibrium phase transition at a finite temperature into the glassy ordered state. 
This has been particularly well established in canonical SGs which represent a class of dilute 
transition-metal alloys soluted in the noble-metal host. 

In theoretical studies of the ordering of SGs, a statistical model called the Edwards- Anderson 
(EA) model has been widely used [2j. In this model, spins are located on a certain regular 
lattice, say, on a three-dimensional (3D) simple cubic lattice, and are assumed to interact via 
the random exchange interaction taking either ferromagnetic or antiferromagnetic sign. One 
crucially important ingredient of the EA model is its spin symmetry. In its original version 
as proposed by Edwards and Anderson, the spin was assumed to be the Heisenberg one with 
three spin components. In fact, this is expected to be a good approximation for most of real 
SG materials, particularly for canonical SGs which possess relatively weak magnetic anisotropy. 
Another well-studied SG model is the Ising EA model with only one spin component. Although 
an ideal Ising SG is relative scarce in real SG materials, it has been widely used in numerous 
SG simulations performed in the past because of its simplicity. A naive expectation might be 
that the weak magnetic anisotropy inherent to real magnets eventually causes a crossover from 
the isotropic Heisenberg behavior to the anisotropic Ising behavior even in real SG materials 
[3], so that the Ising EA model corresponding to the strongly anisotropic limit is enough to 
describe the asymptotic behavior of weakly anisotropic SG magnets. Although such a naive 



expectation sounds plausible from the standard view of critical phenomena, the actual situation 
is not necessarily so simple as we shall see below. 

Obviously, it is also crucially important to fully understand the ordering properties of the 
Heisenberg EA model, in particular, whether the Heisenberg EA model behaves similarly or 
differently from the Ising EA model. Indeed, via recent extensive studies, it now becomes 
increasing clear that the 3D Heisenberg EA model exhibits an intriguing ordering behavior 
quite different from that of its Ising counterpart. Among others, a special degree of freedom 
called a "chirality", totally absent in the corresponding Ising model, might play a crucial role 
in the ordering process of the 3D Heisenberg EA model [1], which might be crucially important 
in the ordering of real weakly anisotropic SG materials [H [3]. 

In the present article, we wish to review the recent studies on the ordering properties of the 
two typical models of SGs, i.e., the 3D Ising EA model versus the 3D Heisenberg EA model, 
focusing on their similarities and differences. Indeed, although the recent numerical studies 
have indicated that both the Ising EA model and the Heisenberg EA model exhibit a finite- 
temperature transition, the nature of the ordering is often in sharp contrast with each other. In 
the present article, we wish to focus on the following four issues. 

(a) Is there a thermodynamic transition in zero field ? 

(b) If so, what are the associated critical properties ? 

(c) What is the nature of the ordered state ? Does it exhibit a replica-symmetry breaking (RSB) 
? If yes, what type ? 

(d) Is there a thermodynamic transition in applied magnetic fields ? 

In the following, we shall examine the ordering properties of the 3D Ising SG in section 2, and 
those of the 3D Heisenberg SG in section 3. Section 4 is devoted to summary and discussion. 

2. The Ising spin glass 

In this section, we review the ordering properties of the 3D Ising EA model. The Hamiltonian 
of the 3D Ising EA model is given by, 

^ = - E JijS^Sj-HY^Si , (1) 

<ij> i 

where Si = ±1 is an Ising spin variable located at the i-th site of a 3D simple-cubic lattice, H 
is an external field intensity, and the < ij > sum is taken over all nearest-neighbor pairs on 
the lattice. Periodic boundary conditions are applied in all directions. The nearest-neighbor 
couplings Jij are quenched random variables with zero mean and standard deviation unity. For 
the form of the Jij distribution, either the Gaussian distribution or the binary distribution (it J 
distribution) is widely used. 

The ordering properties of the 3D Ising EA model has been studied quite extensively for 
years as a standard model of SGs since the beginning of the SG research in 1970's. Note that 
the model has a finite-range (nearest-neighbor) interaction so that, in contrast to the mean- 
field-type (or an infinite-ranged) model such as the Sherrington-Kirkpatrick (SK) model [6], the 
correlation effect is taken into account in the model. In spite of its apparent simplicity, however, 
the ordering properties of the 3D Ising SG model turns out to be highly nontrivial and hard to 
be elucidated. So far, rigorous analytical information on this model has been scarce, and most 
of the information we now have for this model is necessarily a numerical one. Nevertheless, after 
more than 30 years of extensive research, we now have several established (or at least almost 
established) knowledge about the ordering properties of the 3D Ising EA model. We summarize 
some of them below, together with several open issues. 



2.1. Ordering in zero field 

In earlier studies, the issue of whether the 3D Ising EA model exhibits a finite-temperature 
transition or not remained controversial. Monte Carlo (MC) simulations on the model made 
by Ogielski and by Bhatt and Young in 1985 changed this situation, giving strong support to 
the occurrence of a finite-temperature SG transition [71[8]. Since then, many numerical works 
have been made, all of which have agreed in that the 3D EA model exhibits a thermodynamic 
SG transition at a finite temperature [HllinilllllSlIiaillllSlIISlIIZlIISlIIS]- In particular, 
the correlation-length ratio and the Binder ratio played an important role in identifying the 
SG transition point. In Fig.l, the correlation-length ratio ii/L, defined by the finite-size SG 
correlation length divided by the linear size of the system L, calculated by Katzgraber et 
al for the ibJ Ising EA model of L < 32 is shown, together with the Binder ratio calculated 
by the same authors [IS]. The correlation- length ratio and the Binder ratio are dimensionless 
quantities so that the data of various L should be scale-invariant and exhibits a crossing behavior 
at the respective SG and chiral-glass transition points. Indeed, clear fan-out of the correlation- 
length ratio and the Binder ratio is observed at /3~^ = T ~ 1.12, indicating the occurrence of a 
finite-temperature SG transition. 




Figure 1. (Color online) The temperature and size dependence of the correlation- length ratio 
(upper figure), and the Binder ratio (lower figure), of the 3D Ising EA model with the binary 
(± J) coupling. (3 is the inverse temperature. Taken from [H. G. Katzgraber, M. Korner and A. 
P. Young, Phys. Rev. B73 (2006) 224432]. 

By contrast, reliable estimates of the associated SG exponents are more difficult. Concerning 
the critical-point-decay exponent (or the anomalous dimension) ry, various simulations have given 
nearly a common value of ~ —0.38 ~ —0.40. By contrast, concerning the correlation-length 
exponent i/, smaller values around u ~ 1.2 — 1.3 reported by earlier simulations [318] have been 
revised by later simulations to significantly larger values. Indeed, most recent MC simulations 
have yielded i/ ~ 2.5 ~ 2.7, together with r? ~ -0.38 ~ -0.40 O US]. The SG susceptibihty 
exponent 7 is given by 7 ~ 6.0 ~ 6.5. As an example, we show in Fig. 2 the finite-size scaling 
plot of the correlation-length ratio of the it J Ising EA model reported by Campbell et al, which 
yielded v = 2.72(8) [ISj. 

Although a good experimental realization of the Ising SG has been rather scarce, at least 
one example, an insulating SG magnet FeMnTiOs, has been thoroughly studied. Indeed, the 




Figure 2. (Color online) Finite-size scaling plot of the correlation-length ratio of the 3D Ising 
EA model with the binary (±J) coupling. The exponent v is taken to be 2.72. Taken from [I. 
A. CampbeU, K. Hukushima, and H. Takayama: Phys. Rev. Lett. 97 (2006) 117202]. 



experimental exponent values measured for FeMnTiOs turns out to be roughly consistent with 
the numerical values of the 3D Ising SG [201 [^. 

Meanwhile, in many of real SG magnets including canonical SGs, magnetic anisotropy 
is much weaker than the isotropic exchange interaction so that they might be described as 
a weakly anisotropic Heisenberg-like SG. It should be stressed that, in contrast to a naive 
expectation mentioned above, the experimentally determined SG critical exponents of such 
weakly anisotropic Heisenberg-like SGs largely deviate from the exponent values of the 3D 
Ising EA model. Indeed, for typical canonical SG materials, consistent experimental estimates 
are now available thanks to careful experimental measurements [211 [22} [23} [2^ [25} [26} [27] . The 
exponents determined by various authors for canonical SGs indeed come close to each other, 
yielding the values /3 ~ 1, 7 ~ 2.2 — 2.3, ~ 1.3 — 1.4 and r] ~ 0.4 — 0.5, which are clearly at 
odd with the exponent values of the 3D Ising EA model. The experimental v and 7 are about 
half and the sign of r] is reversed. 

To elucidate the origin of the discrepancy observed between the critical properties of 
experimental Heisenberg-like SG magnets and of the 3D Ising SG is one of important issues 
left in the SG research. 

2.2. Possible replica- symmetry breaking 

One of hot issues in the SG research has been concerned with the nature of the SG ordered 
state: In particular, whether the SG ordered state spontaneously breaks a replica symmetry or 
not [T]. Two typical views have been common. One is a droplet picture, which claims that the 
SG ordered state is a "disguised ferromagnet" without a spontaneous RSB [28J. The other is a 
hierarchical RSB picture inspired by the exact solution of the mean-field model, which claims 
that the SG ordered state is intrinsically more complex accompanied with a hierarchical or full 
RSB where the phase space is hierarchically organized in the SG ordered state [29j. Hot debate 
has continued over years concerning which view applies to the ordered state of real SG magnets. 
For more quantitative discussion, it is convenient to introduce an "overlap" variable q, which 
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Figure 3. (Color online) Typical patterns of the overlap distribution function P{(i) in the 
thermodynamic limit. In the case (a) expected in a droplet theory, -P((?) consists of just two 
delta-function peaks at (? = -^qEA- In the case of (b) expected in a hierarchical or full RSB 
picture, -P(g) possesses an additional plateau part connecting the two delta-function peaks at 
q = ig^;^. In the case of (c) expected in a one-step RSB picture, P{q) possesses a central 
peak located at (7 = in addition to two delta-function peaks at g = ^qea- The case (d) is a 
combination of (b) and (c). 



is defined by 

^=l^i:sr^sf\ (2) 

i 

where 5^^"^ represents the i-th spin variable of the "replica" a and the summation is taken over 
all N spins of the system. Replicas a and /3 mean here the two independent copies of the system 
with the same realization of quenched randomness. One can then consider the distribution 
function of the overlap variable P{q) 

P{q') = [<5{q-q')>], (3) 

where < • • • > represents a thermal average and [• • •] an average over the quenched disorder 
(configurational average). Some typical forms of P{q) in the thermodynamic limit is illustrated 
in Fig. 3. 

The droplet picture claims that the overlap distribution describing the SG ordered state to 
be a trivial one in the thermodynamic limit consisting of just two delta functions located at 
q = qEA and at g = —qEA- It means that the SG ordered state consists of unique pure state and 
its symmetry counterpart, irrespective of its apparent complexity in real-space spin pattern. 

In the hierarchical RSB picture, by contrast, P{q) exhibits a continuous plateau part spanning 
between the two delta- function peaks at g = -iiqEA- It means that the phase space is divided 
into infinitely many pure states organized in a hierarchical manner, each of which is separated 
by infinitely high free-energy barrier. 




Figure 4. The overlap distribution function P{q) in the spin-glass ordered state of the 3D Ising 
EA model with the Gaussian coupling. The temperature is T = 0.7 where Tg was estimated as 
Tg = 0.95 lb 0.04. Only the positive side of q is shown. The sizes are L = 4, 6, 8, 10, 12 and 16. 
Taken from [E. Marinari, G. Parisi and J.J Ruiz-Lorenzo, Phys. Rev. B58 (1998) 14852]. 



Other types of RSB have also been known. One well-known example might be a one-step 
RSB, in which P{q) possesses a central (5-function peak at g = in addition to the self-overlap 
peaks at q = ±qEA- In this case, the phase space is divided into many components, but all of 
them, except for itself and its symmetry partner, are completely dissimilar or orthogonal. It 
is realized, e.g., in the ordered states of the mean-field Potts glass or of the mean-field p-spin 
model [H [301 EI]- In the past, such a one-step RSB has often been discussed in the context of 
the structural-glass problem rather than SG problem. The combination of the full RSB and the 
one-step RSB is also possible in certain models. 

Many numerical studies have focused on the 3D EA model to solve this issue of RSB, 
asking whether this model exhibits a hierarchical RSB or no RSB. The controversy has not 
yet been solved completely. However, most of the recent studies have indicated that the overlap 
distribution P{q) persistently exhibits a nontrivial component suggesting the occurrence of 
a hierarchical-type RSB [IH l32l I33j . As a typical example, we show in Fig. 4 the overlap 
distribution P{q) of the 3D Gaussian Ising EA model of various lattice sizes at a temperature 
T = 0.7 which is well below the estimated SG transition temperature Tg = 0.95 ± 0.04 |11] . 

2. 3. Ordering in finite fields 

In contrast to the situation in zero field, the question of the existence/nonexistence of an 
equilibrium transition under external fields has remained quite controversial for years and not 
yet been completely settled. In fact, the question is closely related to the one discussed in the 
previous subsection, i.e., that of an RSB. For example, if the SG ordered state of the 3D Ising 
EA model does not accompany an RSB, the model is likely to exhibit no finite-temperature 
transition under finite fields since applied fields destroy the global Z2 spin-inversion symmetry, 
a symmetry broken at the zero-field transition. If, on the other hand, the SG ordered state of the 
model accompanies an RSB, a finite-temperature transition associated with RSB may become 
possible even in the situation where no global symmetry is left under fields, though it does not 
immediately necessitate the occurrence of an in-field transition. 

Although there is no completely consensus among researchers, most of recent numerical 



studies point to the absence of an equilibrium transition under finite fields [SU [35]. Such 
numerical observations seem consistent with the experimental observation on FeMnTiOs: For 
this compound, Mattsson et al indicated that there was no in-field phase transition [36j, in sharp 
contrast to the zero-field case. Meanwhile, the absence of an equilibrium transition in finite fields 
seem in apparent contrast to the possible occurrence of a hierarchical- type RSB observed in zero 
field. Further study is required to clarify this point. 

3. The Heisenberg spin glass 

In the present section, we deal with the ordering of the 3D Heisenberg EA model. The 
Hamiltonian of the 3D Heisenberg EA model is given by, 

n = -J2 J^A■s,-HY,s! , (4) 

<ij> i 

where Si = {Sf,Sf,Sf) is now a three-component unit vector located at the i-th site of a 3D 
simple-cubic lattice, H is an external field intensity applied along the z-direction, and the < ij > 
sum is taken over all nearest-neighbor pairs. Periodic boundary conditions are applied in all 
directions. The nearest-neighbor couplings Jij are quenched random variables with zero mean 
and standard deviation unity. As in the Ising case, either the Gaussian distribution or the binary 
distribution (it J distribution) is widely used for the form of the Jij distribution. 

3.1. Ordering in zero field 

The ordering properties of the 3D Heisenberg EA model have long been studied [371 EH l39| 

iiQiiiiiiMiiaiMiiiaiMiiiziiisiiigimiia^ 

are hotly debated even now. Earlier numerical simulations on the model suggested in common 
that the model exhibited only a zero-temperature transition [371 133 ESI HOI HI] , in apparent 
contrast to experiments. Common attitude in the community at that time was to invoke the 
weak magnetic anisotropy inherent to real materials to explain this apparent discrepancy with 
experiments, assuming that the weak anisotropy caused a rapid crossover from the Tg = 
Heisenberg behavior to the Tg > Ising behavior [3]. In fact, however, the situation was not 
quite satisfactory as already discussed above, in view of the fact that the experimental exponent 
values measured for canonical SGs are actually far from the Ising SG values, and no clear sign 
of Heisenberg-to-Ising crossover has been observed in experiments. 

In 1992, the present author suggested that the 3D Heisenberg EA model might exhibit a finite- 
temperature transition in its chiral sector 02] • Chirality is a multispin variable representing the 
sense or the handedness of the noncollinear or noncoplanar structures induced by frustration, 
i.e., whether the frustration- induced noncollinear or noncoplanar spin structure is right- or left- 
handed. The scalar chirality x is defined by the product of three neighboring spins by 

X = Si - Sj-x Sk- (5) 

It has subsequently been suggested that, in the ordering of the 3D Heisenberg SG, the chirality 
was "decoupled" from the spin, the chiral-glass order taking place at a temperature higher than 
the SG order, Tcg > Tsg [IHl HH [571 EH [62] . Based on such a spin-chirality decoupling picture of 
the 3D isotropic Heisenberg SG, a chirality scenario of experimental SG transition was proposed 
by the author [321 [M]: According to this scenario, the chirality is a hidden order parameter of 
real SG transition. 

The chirality scenario consists of the two parts, i.e., the "spin-chirality decoupling" part for 
a completely isotropic system and the "spin-chirality recoupling" part for a weakly anisotropic 
system [1]. The first part, the spin-chirality decoupling, is a key ingredient of the scenario. It 
claims that the fully isotropic 3D Heisenberg SG exhibits a peculiar two-step ordering process. 



in which the systems exhibits, with decreasing the temperature, first the glass ordering of the 
chirality at a finite temperature T = Tcg spontaneously breaking only a discrete Z2 symmetry 
with preserving the continuous 50(3) symmetry, and at a lower temperature T = Tsg < Tcg 
exhibits the glass ordering of the spin itself fully breaking both the Z2 and 50(3) symmetries. 
The higher transition at T = Tcg associated with the discrete Z2 symmetry breaking is called 
the "chiral-glass transition" , while the intermediate phase between Tcg Sind Tsg where only the 
Z2 symmetry is broken with preserving the continuous 50(3) symmetry is called the "chiral-glass 
state" . 

The second part of the chirality scenario concerns the role of the weak anisotropy which 
inevitably exists in real SG magnets to certain extent. The chirality scenario claims that the 
weak random magnetic anisotropy, which reduces the Hamiltonian symmetry from Z2 x 50(3) 
to only Z2, weakly "mixes" the chirality to the spin sector, and the chiral-glass transition hidden 
in the chiral sector in fully isotropic system is 'revealed' in the spin sector in weakly anisotropic 
system. In this scenario, the chiral-glass transition of the fully isotropic system, not the SG 
transition of the isotropic system, dictates the SG of real weakly anisotropic SGs. 

The chirality scenario is capable of explaining several long-standing puzzles concerning the 
experimental SG transition in a natural way, e.g., the origin of the non-Ising critical exponents 
experimentally observed in canonical SGs, and remains to be an attractive hypothesis in 
consistently explaining various experimental observations for canonical SGs [4J. In recent 
numerical studies of the 3D Heisenberg EA model, although consensus now seems to appear 
that the 3D Heisenberg SG indeed exhibits a finite-temperature transition of some sort 
[571 EU [581 |59l ESI i63j , the nature of the transition, especially whether the model really exhibits 
the spin-chirality decoupling, is stih under hot debate [IZllMlliaEaEallHllSSllMllMlISS]. 

Thus, the issue of whether the spin-chirality decoupling really occurs in the 3D Heisenberg 
SG has remained controversial in spite of the potential importance in understanding the nature 
of experimental SG ordering. While several numerical results in favor of the occurrence of 
the spin-chirality decoupling were reported [l5l 061 |50l EU [551 EZl IMl [63], a simultaneous 
spin and chirality transition without the spin-chirality decoupling was claimed in other works 
|lZl[iaiia[52l[531[Ml[5a[58l[59l[63]. The recent simulation of Ref.t62j, however, has provided 
a fairly strong numerical support for the occurrence of the spin-chirality decoupling. This 
calculation indicates that the SG transition point Tsg is located about 10% ~ 15% below the 
chiral-glass transition point Tcg^ which we shall briefly describe below. 




Figure 5. (Color online) The temperature and size dependence of the correlation-length ratio 
for the chirality (a), and for the spin (b). The arrow indicates the bulk chiral-glass transition 
point. 
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Figure 6. (Color online) The temperature and size dependence of the Binder ratio for the 
chirality (a), and for the spin (b). The arrow indicates the bulk chiral-glass transition point. 




Figure 7. (Color online) The (inverse) size dependence of the crossing temperatures of E.cg/L 
and S^sg/L, the dip temperature Tdip and the crossing temperature Tcross of gcc- The inset 
exhibits a wider range. 



In Fig. 5, we show the recent MC data of the correlation- length ratios for the chirality ^cg/L 
(a), and for the spin ^sg/L (b), of the isotropic 3D Heisenberg EA model with random Gaussian 
coupling [62]. The system contains total N = spins with L ranging from 6 to 32. For further 
details, refer to Ref . [62] . As can be seen from the figure, while the chiral £,cg/L curves cross at 
temperatures which are only weakly L-dependent, the spin ^sg/L curves cross at progressively 
lower temperatures as L increases. 

As an other indicator of the transition, we show in Fig. 6 the Binder ratios for the chirality 
(a), and for the spin (b). The Binder ratios are also dimensionless, and are expected to exhibit a 
scale-invariant behavior at the respective chiral-glass and SG transition points. As can be seen 
from the figure, the chiral Binder ratio gcG exhibits a negative dip which deepens with increasing 
L. The data of different L cross on the negative side of gcG unlike the correlation-length ratio. 
These features indicate a finite-temperature transition in the chiral sector. 

To estimate the bulk chiral-glass and SG transition temperatures quantitatively, we plot 
in Fig. 7 the crossing temperature of S,cg/L and ^sg/L for pairs of successive L values versus 
1/Lav, where Lav is a mean of the two sizes, together with the dip temperature T^ip and the 
crossing temperature Tcross of the chiral Binder ratio gcG- The data show ci near- linear 1 / Ld^- 
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Figure 8. (Color online) Finite-size-scaling plot of the chiral-glass correlation-length ratio 
^cg/L^ where the correction-to-scaling effect is taken into account. The chiral-glass transition 
temperature and the leading correction-to-scaling exponents are Tqg = 0.143 and uj = 0.27. 
The best fit is obtained with vcg = 1.5. 



dependence. The chiral crossing temperatures of ^,cg/L and of gcG exhibits a weaker size 
dependence than the spin crossing temperature, and are extrapolated to Tcg = 0.143 ± 0.004. 
The spin crossing temperature exhibits a stronger size dependence, which is extrapolated to 
TsG = 0.125 ± 0.006. Hence, Tsg is lower than Tqg by about 10% ~ 15%. 

Various chiral-glass exponents have been estimated via the standard finite-size scaling 
analysis. As an example, we show in Fig. 8 the finite-size scaling plot of the chiral correlation- 
length ratio ^cg/L in which the effect of the leading correction-to-scaling term has been taken 
into account. The exponents determined in this way from various physical quantities yield 
Vcg = 1-5 ± 0.2 and rjcG = 0.6 it 0.2, etc., which differ significantly from the standard 3D Ising 
SG values, i/ ~ 2.5 ~ 2.7 and r] ~ -0.38 ~ -0.40 [ISl HSj. The results indicate that the chiral- 
glass transition belongs to a universality class different from that of the 3D Ising SG. Possible 
long-range and/or many-body nature of the chirality-chirality interaction might be the cause of 
this difference. 

Concerning the critical properties associated with the SG transition, which is likely to be 
located below the chiral-glass transition, the available information is still rather limited. The 
exponent rise was reported to be -qsc ~ —0.3 in Ref. . while the exponent usG to be vsG — 1-5 
in Ref. [63]. 

3.2. Possible replica- symmetry breaking 

In this subsection, we discuss the issue of the nature of the chiral-glass ordered state, i.e., 
whether the chiral-glass state exhibits an RSB, and if so, what type. In this connection, it 
should be noticed that the form of the chiral Binder ratio gcG shown in Fig. 6, which exhibits 
a prominent negative dip, is quite peculiar. In fact, this form of the Binder ratio resembles the 
one of the system exhibiting a one-step RSB [301 [3T]. 

In Fig. 9, we show the chiral-overlap distribution P{qx) in the chiral-glass phase calculated in 
Ref. |57j for the 3D Heisenberg SG model with the binary coupling. Here, the chiral overlap q^. is 

defined by q^, = j^J^if^xl^Jxi^^ where (1) and (2) indicate two copies (replicas) of the system. 
The local chirality variable Xifj. is defined at the ith. site and in the /x-direction (// = x,y, z) 
by Xifi = Si-e^ ■ {Si X Si+e^) where e^ is a unit lattice vector along the //-direction. The 
calculated P{qx) exhibits besides symmetric side peaks located at q^ = i'?^^ corresponding 




Figure 9. The overlap distribution function for the chirality of the 3D Heisenberg spin glass 
with the ibJ binary coupling. The temperature is T = 0.15 below the chiral-glass transition 
temperature of this model T = Tqg — 0.19. [From K. Hukushima and H. Kawamura, Phys. 
Rev. B72 (2005), 144416.] 



to the long-range chiral-glass order, which grow with increasing L, it also exhibits a prominent 
central peak at = 0, which also grows with increasing L. The existence of such a pronounced 
central peak is a characteristic feature of the system exhibiting a one-step- like RSB, never seen 
in the Ising SG. The data strongly suggest that the chiral-glass ordered state exhibits a one- 
step-like RSB [l6l [571 E2], in quite contrast to the one of the 3D Ising SG. Indeed, recent 
off-equilibrium experiments of the fluctuation-dissipation ratio of the Heisenberg SG magnet, 
thiospinel CdCr1.7Ino.3S4, seems consistent with such an expectation [ll[6l]. Recent MC also 
indicates that the chiral-glass ordered state is non-self-averaging \57\ 162). 

3. 3. Ordering in finite fields 

In this subsection, we discuss the ordering of the 3D Heisenberg SG under finite fields [50] . 
When one applies a field to the fully isotropic Heisenberg SG, the symmetry of the Hamiltonian 
reduces from Z2 x 5*0(3) in zero field to Z2 x SO (2) under fields, where Z2 refers to the chiral 
degeneracy associated with a spin-refiection operation (solely in spin space, not in real space) 
with respect to an arbitrary plane in spin space including the magnetic-field axis, while S0{2) 
refers to the continuous degeneracy associated with a spin-rotation operation (in spin space, not 
in real space) around the magnetic-field axis in spin space. 

Since the Z2 chiral symmetry characterized by the sign of the scalar chirality remains under 
magnetic fields, the chiral-glass transition is expected to persist under magnetic fields. Of 
course, applied fields change the symmetry, but lower only the continuous part from SO (3) 
to SO (2). Since the continuous part has already been decoupled from the discrete Z2 part, 
a natural expectation here would be that the chiral-glass transition persists even under fields 
essentially of the same type as the zero-field one. In particular, the chiral-glass transition line 
under fields should be a regular function of the filed intensity H. Since there is a trivial H —H 
symmetry, the chiral-glass transition temperature under fields should behave for weak fields as 
Tcg{H) PS Tcg(O) — cH^ • • • (c is a constant). In fact, this yields a transition line resembling 
the so-called GT line of the mean-field model [05], \TcGiO) - Tcg{H)\ oc F^/^ although the 
origin of the exponent 1/2 is entirely different: Here, 1/2 is just of regular origin, whereas the 
exponent 1/2 in the mean-field model cannot be regarded as of regular origin. The expected 




Figure 10. Typical phase diagram of the 3D isotropic Heisenberg SG in the magnetic field 
(H) - temperature (T) plane, expected from the chirality scenario. 'P', 'SG' and 'CG' represent 
the paramagnetic, spin-glass and chiral-glass phases, respectively. 



phase diagram expected from the chirality scenario is sketched in Fig. 10. 

The ordering associated with the continuous part should also occur under finite fields. The 
in-field transition line associated with the continuous 50(2) symmetry breaking should be a 
continuation of the 50(3) breaking SG transition in zero field. Since the broken symmetry 
is different in zero and finite fields, i.e., 50(3) versus SO{2), the SG transition line at low 
fields should exhibit a singular form, H oc |r5G'(0) — TsGiH)\'^^'^ . The exponent (p is not yet 
precisely determined, but may roughly be estimated as i;^ = I^SG + ISG ~ 4. This 50(2) breaking 
transition line is also sketched in Fig. 10. 

In-field ordering properties of the isotropic 3D Heisenberg SG were studied by Imagawa 
and the author also by means of MC simulations, which successfully demonstrates some of the 
feature of the phase diagram mentioned above |501 155j. The obtained magnetic phase diagram 
is also consistent with the experimental one for canonical SGs O [66| I67j . In Fig.ll, we show 
the magnetic phase diagram of the 3D isotropic Heisenberg EA model with binary coupling as 
determined by MC simulation [50]. (Note that the SG transition line was not examined in this 
work). Another interesting observation from MC is that the SG ordered state turns out to be 
quite robust against applied magnetic fields as can be seen from Fig.ll [501 EH]- It appears to 
be stable up to fields as large as 25kBTsG{H = 0). This might be understandable if one notices 
that the coupling between the chirality and magnetic fields might be rather weak, since magnetic 
fields couple directly to the spin via the Zeeman term, only indirectly to the chirality. 

4. Summary and discussion 

In this article, I reviewed the ordering properties of the two representative models of SGs, the 
3D Ising EA model and the 3D Heisenberg EA model. Particular attention has been paid to i) 
whether the model exhibits a finite-temperature transition in zero field or not, and if it does. 
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Figure 11. Magnetic phase diagram of the 3D isotropic Heisenberg SG in the magnetic field 
(H) - temperature (T) plane. Note that the energy scales of the magnetic field H and the 
temperature T are mutually different in (a), while they are taken to be common in (b). Taken 
from [D. Imagawa and H. Kawamura, J. Phys. Soc. Jpn. 71 (2002) 127]. 

what are the critical properties of the transition; (ii) what is the nature of the ordered state, 
whether the model exhibits an RSB, and if yes, what type; (iii) whether the model exhibits a 
finite-temperature transition in finite fields. 

Concerning the zero- field ordering, both models, i.e., the 3D Ising and Heisenberg EA models, 
exhibit a finite-temperature transition. A single SG transition occurs in the Ising EA model, 
whereas successive chiral-glass and SG transitions are likely to occur in the Heisenberg EA model, 
though there still remains controversy about the occurrence of such spin-chirality decoupling. 
Concerning the associated critical properties, recent numerical estimates yield use = 2.5 ~ 2.7 
and r]sG = —0.38 ~ —0.40 for the 3D Ising SG, while, for the chiral-glass critical properties 
of the 3D Heisenberg SG, recent numerical estimates yield vcG ~ 1-5 and rjcc ~ 0.6. For the 
SG exponents of the 3D Heisenberg SG at its SG transition, which is likely to be located below 
the chiral-glass transition, recent numerical calculation suggests use ~ 1-5 and rjsc ~ —0.3. 
One sees from these results that the universality classes associated with the spin ordering of the 
Ising SG, the chiral and the spin orderings of the isotropic Heisenberg SG are all different from 
each other. Numerically observed critical properties of the 3D Ising SG are roughly consistent 
with the experimental values of the Ising-like SG magnet FeTiOs. Numerically observed chiral- 
glass critical properties of the 3D Heisenberg SG are close to those experimentally observed 
for canonical SGs. This coincidence supports the chirality scenario of Refs.|42j which claims 
that the SG critical properties of real weakly anisotropic Heisenberg-like SG magnets such as 
canonical SGs should be equal to the chiral-glass {not SG !) critical properties of the fully 
isotropic Heisenberg SG. 

Concerning the issue of the RSB in the ordered state, numerical simulations on the 3D Ising 
SG point to the occurrence of a hierarchical- type RSB like the one in the infinite-range SK 
model, although the situation has remained somewhat controversial. For the 3D Heisenberg 
SG, recent simulations suggest that the ordered state might exhibit a one-step-like RSB, the 
one quite different from that of the Ising SG or of the mean-field SK model. According to the 
chirality scenario, such a one-step-like RSB should be realized in real weakly anisotropic SG 
magnets. 

Concerning the in-field ordering properties, recent numerical simulations suggest that the 3D 



Ising EA model does not exhibit a finite-temperature transition in magnetic fields, which seems 
consistent with the experiment on the Ising-like SG magnet FeMnTiOs. By contrast, the 3D 
Hciscnbcrg EA model exhibits a finitc-tcmpcraturc transition in magnetic fields associated with 
its chiral degrees of freedom, consistently with experimental observations for canonical SGs. 

Thus, the ordering properties of the 3D Ising SG and the Heisenberg SG are quite different 
from each other, each of which might provide a useful reference in attacking various problems 
related to SGs, not only of random and/or frustrated magnets but also of many systems including 
glassy materials, complex systems and even information problem. 
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